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b. 

We present the results of a search for the three neutral charm decays, D° — » fx e^, D° — » /i + fi , 
and D° — * e + e~ . This study was based on data collected in Experiment 789 at the Fermi National 
Accelerator Laboratory using 800 GeV/c proton-Au and proton-Be interactions. No evidence is 
found for any of the decays. Upper limits on the branching ratios, at the 90% confidence level, 
of 1.56 x 10~ 5 for D° -> ^+ pT , 8.19 x 10~ 6 for D° -> e+e" and 1.72 x 10" 5 for D° -> /j^e* are 
obtained. 

PACS numbers: 13.30.Fc, 13.25.Ft, 13.85.-t, 14.40.Lb 



I. INTRODUCTION 

In the Standard Model, the flavor-changing neutral-current decays D° — > e + e~ and D° — * n + fi~ are forbidden at 
tree leveljj At the one-loop level, the decays are GIM- and helicity-suppressed. Thus the branching ratios are expected 



L In this paper, the symbol D denotes both D and D mesons. 



to be very small. Lepton-family-number violating decays such as D° — > /x ± e =F are strictly forbidden. However, 
extensions of the Standard Model allow for both flavor-changing neutral currents and lepton-family-number violation, 
so detection of such dilepton decays could be taken as evidence for new physics (l]-|| . Previous experiments || ||] 
have quoted limits of order 10 -4 to 1CP 6 for several such decays. In this paper we describe an experiment that places 
limits between a few times 10 -5 and 10 -6 on the branching ratios for the decays D° — > /i + /i~, D° — > e + e~, and 
D° -» ^e* 

II. THE E789 SPECTROMETER 

Experiment 789 was carried out in the Meson East beam line at the Fermi National Accelerator Laboratory, where 
a beam of 800 GeV/c protons was delivered to a fixed target of either gold or beryllium. 

The spectrometer, shown in Figure |l|, was optimized for two-body final states with pair rapidity near zero in 
the center-of-mass system. Its main components were a silicon-strip vertex detector (SSD) just after the target, a 
copper beam dump, two dipolc bending magnets (SM12 and SM3), three stations of drift chambers and hodoscopes, 
a sampling calorimeter, and a muon-identification station located at the end of the spectrometer. The ring-imaging 
Cherenkov detector was not used in this analysis. 

A. Target 

The target apparatus was installed in the beam vacuum. Table p] gives the dimensions of the targets used for this 
analysis. The targets were much wider than the beam in the x (horizontal) dimension but were narrow in the y 
(vertical) dimension. The target centers were located at z = —331.85 cm. Here, the z axis is the direction of the 
incident proton beam, and the origin of the right-handed coordinate system is centered at the upstream end of the 
SM12 yoke. 

B. Beam Monitors 

Beam intensity was measured using both an ion chamber and a secondary-emission monitor, SEM, located upstream 
of the target. The fraction of beam striking the target was determined using an interaction monitor, AMON, which 
was a scintillation-counter telescope perpendicular to the beam and viewed the target through a hole in the shielding 
cave. The targeting fraction varied from 30% to 40% depending on the running conditions (see J7| for details). 

C. Silicon Vertex Detector 

The SSD was located just downstream of the target and consisted of two arms, each containing eight planes of 
detectors (see Figure ||). Each 5-cm by 5-cm plane was a 300-^m-thick silicon-strip detector with 50- ^m strip pitch. 
The planes were arranged in two arms to cover vertical angles from 20 to 60 mr above and below the beam. Each 
plane had one of three orientations, Y, U, or V, with rotations about the z axis of 0°, +5°, or —5°, respectively. The 
sequence of orientations in each arm was YUYVYUYV, proceeding downstream. 

A total of 8,544 strips were instrumented with amplifiers S, discriminators 0, and latches [Gil having «30-ns 
effective time resolution. To minimize secondary interactions, thermal fluctuations, and radiation-induced detector 
degradation, the SSD containment volume was temperature-controlled with a 10°C helium fill. Table |l| gives the 
configuration of the SSD planes. One 1-mm-thick scintillator, of dimensions 5 cm x 5 cm, was placed at the downstream 
end of each arm and was used for triggering. 

D. Beam Dump 

A water-cooled copper beam dump was located inside the SM12 magnet (see Figure 0). It prevented noninteracting 
primary protons and secondary particles of low transverse momentum from entering the downstream spectrometer. 
The tapered dump and the baffles on the inside walls of SM12 defined the spectrometer aperture. 



E. Spectrometer Magnets 

The two dipole magnets, SM12 and SM3, served to focus charged particles of momenta within a desired range 
onto the downstream detectors. The magnetic fields of both magnets were carefully mapped with the Fermilab 
Ziptrack [pTj, and the resulting profiles of the field were used in the data analysis. 

SM12 was a 1200-ton, 14.5-m-long, open-aperture dipole magnet. The horizontal aperture was tapered to provide 
a gradually decreasing magnetic field. At an operating current of 900 A, SM12 provided a vertical transverse impulse 
(p t ) of 1.6 GeV/c, optimal for studying a D° decaying into two charged particles. 

The second bending magnet, SM3, was located between tracking stations 1 and 2. It was a 3.4-m-long, open-aperture 
magnet and provided a 0.91 GeV/c vertical pt impulse. It deflected charged particles in the opposite direction as 
SM12, focusing them onto the subsequent detectors. In combination with the drift chambers SM3 provided momentum 
analysis for charged particles. 

F. Tracking Stations 

Three drift-chamber tracking stations were used to determine charged-particle trajectories through the spectrometer. 
Each station consisted of three pairs of chambers, with the chambers in each pair offset by half a drift cell to resolve 
the "left-right" tracking ambiguity. Each pair was oriented in one of three views, Y, U, and V, at 0°, +14°, and —14° 
with respect to the y axis. The drift gas was a 50/50 mixture of argon and ethane, with a 0.7% admixture of ethyl 
alcohol. Each sense wire was connected to its own time-to-digital converter to measure the drift time of the ionization 
electrons. At operating voltages around 2000 V, the drift velocities averaged about 50 /im/ns. Hodoscope planes at 
each tracking station provided fast coarse tracking information used in the trigger. Stations 1 and 3 had both X and 
Y hodoscope planes (designated HY1, HX1, HY3, HX3), while station 2 had only a Y plane (HY2). Each hodoscope 
plane consisted of two half-planes of scintillation counters, whose light was collected using lucite light guides glued to 
Hamamatsu R329 photomultiplier tubes. 

G. Calorimeter 

Sampling calorimeters were used to identify electrons and hadrons (see Figure [|). The electromagnetic section 
consisted of four lead/scintillator layers, El, E2, E3, and E4, with thickness of 2, 5, 5, and 6 radiation lengths, 
respectively. The total thickness of the electromagnetic section was 0.81 interaction length. Each layer had separate 
left (x > 0) and right (x < 0) sections which were divided into twelve modules in y. Each of the resulting 96 modules 
was read out individually, with the analog signal from the scintillator converted to a digital signal using an 8-bit 
quadratic ADC @. 

The hadronic section consisted of two iron/scintillator layers, HI and H2, of 2.14 and 5.84 interaction lengths 
respectively. The left and right segments of each were divided into thirteen modules in y, giving a total of 52 modules. 
For details regarding the calibration of the calorimeters see |13|] . 

H. Muon Station 

The Muon Station, located at the downstream end of the spectrometer, contained three planes of proportional- 
tube arrays and two planes of hodoscopes, interspersed with shielding. The proportional tubes, used in the Trigger 
Processor and for off-line muon identification, consisted of two planes of horizontal cells, PTY1 and PTY2, for y 
determination and one of vertical cells, PTX, measuring the x coordinate. Each plane was made of a series of two- 
layer aluminum extrusions, each containing fifteen 2.54 cm x 2.54 cm cells. The layers were offset by half a cell width 
from each other. The cells were read out with latches, so no timing information was recorded. The gas mixture used 
was the same as in the drift chambers. 

The muon hodoscopes were used both for triggering and for particle identification. There were two planes, HY4 and 
HX4, providing y and x information respectively. The calorimeter and additional zinc, lead, and concrete shielding 
comprised 16 interaction lengths of material between station 3 and the muon station. In addition, concrete absorbers 
interspersed among the muon detectors added approximately 5 more interaction lengths of shielding. 



III. DATA ACQUISITION 



The data-acquisition system um was based on the Nevis Laboratories Data Transport System pa] . Event informa- 
tion from the front-end crates was buffered into multiport memory modules and supplied to the Trigger Processor |lq | 



(described in Section IV C) before readout to the VME-based archiving system, which recorded data on four Exabyte 
8200 tape drives. The system was capable of streaming approximately 1 MB/s to tape. Once per spill, scalers were 
read out to record trigger rates and beam-intensity information with and without system deadtime. 

IV. TRIGGER 

E789 utilized a three- level trigger system. Level- 1 triggers based on hodoscope information were OR'ed together to 
form the Trigger Fan In (TFI) signal. Events satisfying TFI were latched and fed to the "DC Logic" system, in which 
further logical requirements were imposed. In the DC Logic, information from slower spectrometer components, such 
as the calorimeter, could be included. Events satisfying the DC Logic requirements produced the Trigger Generator 
Output (TGO) signal, which vetoed the fast latch reset, preserving hit information for readout to the Trigger Processor. 
Finally, the Trigger Processor examined hit patterns in the wire chambers, hodoscopes, calorimeter, muon detectors, 
and silicon detectors to enhance the fraction of events in the desired decay channels that contained decay vertices 
downstream of the target. Events satisfying all three levels of trigger were written to tape. In addition, at each trigger 
level, some events were prescaled and forced through to the next level. 

A. TFI 

The TFI had three main components, designed to trigger on pairs of charged particles from the target. The main 
pair trigger, |M, required at least two triple hodoscope coincidences in HY1, HY2, and HY3, each corresponding to a 
different charged-particle trajectory from the target. As shown in Figure pi these "Trigger Matrix" coincidences were 
implemented as a look-up table, using fast ECL RAM, that provided four independent output signals corresponding 
to hodoscope roads to the left or right of the z axis and passing above or below the beam dump. At least two of these 
four outputs needed to fire to satisfy %M. 

To allow sufficient redundancy for hodoscope and trigger efficiencies to be determined off-line, additional triggers 
were implemented using majority logic on combinations of the six hodoscope planes HXI, HY2, HY3, HX3, HY4, 
and HX4. These were designated jfiLR and jLR. The notation j/iLR represents a logical AND of j fiL with \[i.R. 
The component j/j,L required that at least n of HXI, HY2, HX4 and HY4 had hits to the left of the z axis. Likewise 
\[iR required at least n from the same group to have hits on the right side. The trigger jLR imposed a similar 
requirement except that planes HXI, HY2, HX3, and HY3 were used. In the run with SMI2 current set to 1000 A, 
n was set to 3. It was changed to 4 in the 900 A run. 

B. DC Logic 

The DC Logic, the second-level trigger, incorporated information from slower detectors whose use at the TFI stage 
would have imposed excessive deadtime. To reject tracks missing the SSD planes, signals from the scintillators behind 
the SSD arms, SU and SD, were required at this stage. Veto signals (NX1 and NX3), formed by counting the 
number of hit counters in the hodoscopes HXI and HX3, were used to veto high-multiplicity events. The DC Logic 
also included particle-identification components based on the calorimeter and the muon station. The various logic 
combinations were OR'ed together to form the TGO signal. 

The DC Logic event-identification requirements caused differences in acceptance for various types of events. The 
dimuon TGO component (/i + /i~) required that 2HXA and 2HYA be satisfied, that is, that two counters in each 
of the muon hodoscopes HX4 and HY4 register a hit. The calorimeter provided a sum of analog signals from the 
dynodes of the E2 and E3 photomultiplier tubes which was sensitive to electrons, and another sum based on signals 
from HI, H2, El and E4 for hadrons. Each sum was discriminated with a low and a high threshold separately. The 
low threshold was set for detecting single particles and the high threshold for two-particle events. The discriminated 
signals e, E, h, and H represented the low and the high thresholds for the electromagnetic and hadronic ener gy s ums 



respectively. Table |Hj lists the various logical combinations that together formed TGO. As shown in Table 



[V, the 



DC Logic requirements reduced the trigger rate significantly relative to TFI. 



C. Trigger Processor 

If an event satisfied one of the DC Logic triggers, the Trigger Processor then searched for tracks from the target 
using the drift-chamber information. Only wire positions were used at this stage. Track hypotheses were formed from 
hits in the Y drift chambers in stations 1, 2, and 3, masked by hodoscope and calorimeter or proportional-tube hits. 
These potential target tracks were then projected to the SSD, and used to identify SSD hits in the y-z view for SSD 
trackfinding. SSD tracks formed from the masked hits were subjected to the requirement that the impact-parameter 
be more than 51 ^m from the center of the target, designed for finding tracks coming from D° decays occurring 
downstream of the target. The chosen tracks were then combined in up-down pairs to form vertices. A cut of 0.10 cm 
was made on the location of the vertex in z to further increase the likelihood that the event contained a downstream 
decay. 

The Trigger Processor reduced the trigger rate by about an order of magnitude below the TGO rate. The Triggers 
After Processor (TAP) was dominated by the dihadron trigger. During the "Dedicated Dilepton" running period, 
the dihadron trig ger was prescaled by a factor of 32 and the proton intensity was increased to enhance the dilepton 



sensitivity. Table IV gives the average rates per spill for protons on target, TFI, TGO, and TAP. 



V. DATA SETS 

Data was taken initially with SM12 set at 1000 A and two different target materials for studying the nuclear 
dependence of proton-induced charm production [|17| . Subsequently data was collected at 900 A because of improved 
acceptance for detecting charm decay at this setting. Three data sets are included in this analysis: lOOOA-Au, 
900A-Au, and 900A-Be. Each set was processed separately and each yielded an independent normalization signal 
in the D° — > Ktt mode. (The latter part of the 900A-Au sample, for which the dihadron trigger was prescaled as 
just described, is referred to as the Dedicated-Dilepton run, but it shared a common normalization signal with the 
rest of the 900A-Au sample.) Table |v| gives the total number of protons on target for each sample, the number of 
AM ON ■ SB countsn (number of live-time-corrected counts in the targeting monitor), and the number of triggers 
recorded on tape. 

VI. NORMALIZATION APPROACH 

To determine the branching ratios for D° — > e + e~, D° — ► /i e^, and D° — > n + fi~, we need the total number of 
D°s produced, as measured by the decay mode D° — > Ktt, as well as the detection efficiency for each decay mode: 

B(D° - l+l~) = ^- x ^B(D° -» Ktt) (1) 

E£X J± xB (J)0, fo) . (2) 

Here Ni+i- is the number of D° — > l + l~ events seen, Nk-k is the number of D° — > Ktt events, £;+;- is the efficiency for 
observing a D° — » l + l~ event, £ktt is the efficiency for observing a D° — ► Ktt event, B(D° — » Ktt) = (3.85 ± 0.09)% 
is the branching ratio for D° — > Ktt decay [18], and e = ck-k/^i+i- is the relative efficiency. Using D° — ► Ktt as a 
normalization mode allowed partial cancellation of many common correction factors in the efficiency ratio. 

VII. EVENT RECONSTRUCTION 

With a total nuclear inelastic cross section per nucleon of 17 mb for beryllium and an inclusive D° production 
cross section of rj 40 jih at ^/s = 39 GeV JLJj], a branching ratio B(D° — * Ktt) of w 4% implies a search for one 
normalization event per ten thousand interactions. In this experiment, with an average D decay distance of rj 3.4 mm 



2 AMON is proportional to the number of interactions in the target and is described in Section [IB. SB is true when the 



system is able to accept data. J^/q'n is thus a measure of the live time of the data acquisition and was typically about 50%. 



(corresponding to an average D° momentum of 56 GeV/c), precise reconstruction of the decay vertex in the SSD is 
a powerful tool to separate D° decays from background processes that occur in the target. The SSD allowed precise 
reconstruction of the decay distance and the impact parameter for each track. The impact parameter could be used 
to eliminate tracks from the target and thus reduce the background significantly. 

A. Pass One 

In the first pass of data processing particle trajectories in both the downstream spectrometer and the SSD were 
reconstructed from the raw data. 

Downstream track reconstruction began by finding hit clusters in the drift chambers. Track segments formed in 
stations 2 and 3 were projected to station 1 for confirmation that the track came from the target and not the beam 
dump. An 18-planc (3 stations with 6 chambers each) least-squares fit was performed. 

The track momentum was determined from the bend angle through SM3. The momentum resolution of the spec- 
trometer was found to be a/p = 1.58 x 10~ 4 p, where p is the momentum of a track and a is the statistical uncertainty. 
The track was then traced back iteratively through SM12, through the SSD, to the target. On the first iteration 
the track was traced from SM3 to the z location of the target center. Candidate tracks falling within an aperture of 
±12.7 cm from the target center in x and y were kept. In subsequent iterations, the track parameters were adjusted 
so that the track traced back to the target center. 

After downstream tracking, all track segments in the SSD were reconstructed. In each of the SSD arms, the four Y 
planes were used first. The preliminary y-z tracks were then employed to define windows in the U and V planes for 
selecting hits that were used to form SSD tracks in the x-z view. Those SSD tracks with enough hits in the Y, U and 
V views were fit to straight lines in three dimensions. The resolution of the impact parameter in y was determined to 
be 34 /im. 

For each event, each opposite-sign pair of downstream tracks was reconstructed as though they decayed from a 
single parent through each of the decay modes D° — > Kir, D° — > /i + /i~, D° — > /i ± e T , and D° — > e + e _ . For at least 
one of these modes, the resulting invariant mass was required to fall within a 500 MeV/c 2 window extending from 
1.65 GeV/c 2 to 2.15 GeV/c 2 . In addition, events were required to have at least one opposite-sign pair of SSD tracks 
that formed a vertex with z location outside a ±2.55 mm window centered at the target. Since the resolution of the 
vertex in z was about 0.7 mm, this requirement rejected a significant fraction of the dihadron events originated from 
the target but still retained about 50% of the D° decays. This pass provided a four-to-one data reduction from the 
raw data. 



B. Pass Two 

In this pass, to further reduce the number of unwanted SSD tracks, the y hits, y-z and x-z angles of the SSD tracks 
were required to be within ±0.2cm, ±1.5 mr and ±2.85 mr respectively from the projections of the downstream tracks 
at the SSD. Figures |6| and [7| show the matching in y and in track angles respectively, before the cuts, for events with 
only one SSD track in either arm. 

The matched SSD tracks, one from each arm, were combined to form pairs. The x /degree of freedom was min- 
imized for each pair by adjusting the vertex location and the track orientation. The downstream tracks were then 
iteratively traced back to the decay vertex as determined by the SSDs to improve the resolution of the track angles. 
The requirement on the invariant mass of the event was tightened to a 200 MeV/c 2 window about the D° mass 
(1.864 GeV/c 2 [Q) for at least one of the modes. This second pass provided another factor of five reduction of the 
data set. 



C. Pass Three 

The alignment of the SSD was refined in this pass. The change in the alignment constants was insignificant. At 
this point, some of the events still contained multiple vertices that resulted either from multiple pairs of downstream 
tracks or from multiple SSD tracks matching a single downstream track. Events were then excluded if more than 
four SSD tracks matched either downstream track or more than ten vertices were reconstructed. To select the proper 
vertex in the surviving events, the quality of each vertex was evaluated using nine parameters. The value of each 
parameter was converted to a probability, with the overall probability taken as the product of the nine probabilities. 
The nine parameters were: 



• % 2 /degree of freedom for reconstructing each SSD track. (2 parameters) 

• x 2 /degree of freedom for the SSD vertex-constrained fit. (1 parameter) 

• y-angle match between each downstream track and its SSD track. (2 parameters) 

• x-angle match between each downstream track and its SSD track. {2 parameters) 

• x 2 /degree of freedom for the position difference between the projection of each downstream track and the SSD 
hit at each SSD plane. (2 parameters) 

Events with only one SSD track in an arm were used to obtain the standard deviations of the distributions for the 
x and y angle matching, as 0.95 mr and 0.25 mr respectively. Only the vertex with the highest overall probability, 
along with the associated SSD tracks and downstream tracks, was employed in the subsequent analysis. 

In the final stage of event selection, a fully reconstructed event consisted of one opposite-sign pair of SSD tracks 
that matched one pair of downstream tracks. The most effective variable to optimize the D° ~ + Kit signal was the 
impact parameter of each SSD track with respect to the target center in y before the vertex-constrained fit. Incorrectly 
reconstructed events often contained at least one track originating in the target that was thus reconstructed with small 
impact parameter. In addition, a cut was made on the lifetime significance, defined as the ratio of t he z location of 
the vertex to the average decay distance of the D° in the laboratory frame, -fftcr. Tables VI, VII and VIII summarize 
the requirements on the impact parameter and lifetime significance for all data sets. 

VIII. PARTICLE IDENTIFICATION 
A. Electron and Hadron Identification 

Electrons and hadrons were identified using the calorimeter. The identification procedure included two requirements: 
first, that energy deposited in the calorimeter match a track, and second, that the profile of the energy deposition in 
the calorimeter be consistent with either a hadron or an electron. 

For each event, the ADC counts of all calorimeter modules were converted to energy. The reconstructed particle 
trajectory was then projected through each layer of the calorimeter. At each layer, the energies deposited in the 
module on the trajectory and in its nearest neighbors were summed. A correction was applied for attenuation of 
scintillation light in the x (readout) direction. If the track had no other track within two modules at each longitudinal 
layer, it was considered isolated and its total energy as well as the "EM fraction" was recorded. The EM fraction 
is the amount of energy deposited in the electromagnetic portion of the calorimeter divided by the total deposited 
energy. 

The energy resolution of the calorimeter was derived from the ratio of the deposited energy in the calorimeter to the 
magnetically- measured momentum of isolated tracks {E/p). R For the hadronic part of the calorimeter, the resolution 
was a/E = -0.018 + 0.91/VE. The energy resolution of the electromagnetic calorimeter was measured to be a jE = 
-0.04 + 0.79/VE. 

For an isolated track, the particle associated with the track was labeled as an electron if the EM fraction was > 0.95 
and the E/p of the track was within 2.58a from the mean of the E/p distribution for that energy bin. A separate 
study, using J/ip — > e + e~ decays from data collected in an adjacent running period, found this cut to be > 96% 
efficient M. For hadrons, the EM fraction was required to be less than 0.7 and the 2.58ct E/p cut was used. 

If two non-muon tracks shared at least one module (which happened quite rarely) , they could often be identified 
using the energy deposition in each section (EM and hadronic) of the calorimeter separately. The procedure for 
identifying the overlapping tracks is described in |13| . 

B. Muon Identification 

Muon identification depended primarily on the muon station. A track was projected to the muon station and each 
detector plane was checked to see if there were hits in momentum-dependent hit windows. The windows, 3ct wide, 



3 The E/p distributions were energy-dependent with various mean and a. These differences were taken into account in the 
analysis as described in J13fl. 



were determined from fits to residual distributions with respect to the projected track. For a track to be categorized 
as a muon candidate, hits matching the projected track were required in both muon hodoscope planes and in at least 
two proportional-tube planes. The hodoscopes had a time resolution better than one 19-ns accelerator-RF bucket. 
Requiring them to have hits dramatically reduced the number of out-of-time tracks. 

In addition to the muon station, the calorimeter was also employed for muon identification. 99.99% of muons under 
100 GeV/c left less than 35% of their energy in the calorimeter. If a track passed the muon hit criteria and had E/p 
greater than 35%, it was tagged as ambiguous but still counted as a potential muon. The most likely mechanism for 
muons to have high E/p was for them to overlap with a non-muon track in the calorimeter. To remove fake dimuon 
events due to a non-muon overlapping with a muon, an isolation requirement was applied. In this case, each track in 
the reconstructed dimuon was required to have unique muon-hodoscope hits and no more than one proportional-tube 
plane could contribute a shared hit. 

IX. ACCEPTANCE AND EFFICIENCY 

It is necessary to calculate the ratio of (acceptance x efficiency) for each dilepton decay to that for the normalization 
decay. This relative acceptance depends on both the kinematics and particle types in each decay mode. While track- 
reconstruction efficiencies cancel in the ratio, trigger and particle-identification efficiencies do not and are determined 
as described below. 

A. Monte Carlo Simulation 

The Monte Carlo (MC) program, incorporating trigger and particle-ID information, was used to generate events 
in each dilepton mode as well as the normalization mode. The MC used the same alignment constants and magnetic 
field maps as the data analysis. Multiple scattering including non-Gaussian tails was also included in the simulation. 
Detector efficiencies were included, and noise hits in the silicon detectors (extracted from data) were added to each 
generated event. 

The production of D° by an 800-GeV proton beam was simulated with a longitudinal fractional-momentum (xf) 
distribution of the form (I — |ie,f|) 6 ' 9 |19|| . The transverse-momentum (pt) distribution was characterized as e~ bPt 
with b — 0.84 (GeV/c) -2 |19| . Each two-body D° decay was generated with a uniform angular distribution in the 
rest frame of the D°. After boosting to the laboratory frame, the decay products (Kir, [ifi, /ze, or ee) were traced 
through the simulated geometry of the spectrometer. Kaon decay was also included in the Monte Carlo. Each event 
that passed the geometric restrictions was required to satisfy the trigger as modeled for the decay of interest. Figure @ 
shows the generated and accepted distributions in p t , xp, and momentum of the D° in the laboratory frame at 900 A. 
With the kinematics of each decay properly modeled, the acceptances were calculated using 40,000 Monte Carlo events 
that passed the geometric cuts for each mode. 

B. Dimuon Efficiency 

Muons were identified primarily using the muon hodoscopes and proportional tubes. The efficiency of each 
proportional-tube plane was determined separately from data and then included in the Monte Carlo simulatio n. Each 



muon selected for the efficiency study satisfied the muon identification requirements as described in Section VIII B . 
This study determined the average efficiencies of the proportional-tube planes to be 95%, 98%, and 95% for Yl, X, 
and Y2 respectively. 

The only component of the dimuon trigger that was not included in the dihadron triggern was the requirement 
of hits in at least two of the four quadrants (upper-left, upper-right, lower-left, and lower-right) in both HX4 and 
HY4. Determining the dimuon efficiency thus required an unbiased study of the muon-hodoscope efficiencies. A muon 
sample was chosen by selecting events that satisfied the calorimeter trigger and included at least one reconstructed 
muon, selected by requiring hits in at least two of the three proportional-tube planes and at least one muon hodoscope 
plane. The window for finding the hodoscope hit was identical to the momentum-dependent window of the nearest 



4 The dihadron-trigger efficiency was used in the calculation of the relative trigger efficiency between the D° — > n + n decay 
and the D° — » Kir decay. 



proportional-tube plane. The requirement of a muon hodoscope plane, whose timing had single-bucket resolution, 
assured that the muon was indeed associated with the triggered event. The efficiency of each hodoscope plane was 
then determined by recording the fraction of events for which the hodoscope that was not used in the muon-selection 
process fired. The efficiencies of HX4 and HY4 were determined to be 95% and 92% respectively, independent of 
muon momentum. These efficiencies were then included in the MC simulation. 

To avoid misidcntification from overlapping tracks, an isolation criterion using the proportional tubes was applied. 
The trigger requirement already isolated the tracks such that the additional proportional-tube isolation criterion 
reduced the efficiency by only 7% while reducing the background significantly. The overall dimuon efficiency was 36% 
at 900 A and 50% at 1000 A. 

C. Dihadron Efficiency 

The efficiency of detecting the D° — ► Kir decay relative to that for the dilepton modes is dominated by the efficiency 
of the H trigger component, which required a significant amount of energy deposited in the calorimeter. 

An unbiased sample of events passing the TFI trigger was employed for studying the H efficiency. Each event in 
this "prescaled-TFI" sample had two hadron tracks with a reconstructed Kit invariant mass in a 500-MeV/c 2 window 
about the D° mass. The momentum range of the selected events was similar to that of the accepted D° — > Kir events. 
The fraction of the prescaled-TFI events that fired the dihadron trigger was then plotted as a function of the total 
momentum in 2 GeV/c bins. Figure O shows the efficiency of the dihadron trigger as a function of momentum and 
the fit thereto by a third-order polynomial. This efficiency curve was then input to the Monte Carlo. The average 
dihadron trigger efficiency for D° — > Kir events that passed the geometric acceptance of the MC was 55% at 900 A 
and 58% at 1000 A. 

The energy deposition of hadrons in each section of the calorimeter was also included in the MC simulation. To 
enhance the certainty of hadron identification, the EM fraction of the MC events was required to be less than 70%. 
This cut accepted over 92% of hadrons. Furthermore, the E/p of the particle was required to fall within ±2.58cr of 
the mean. 

Kaon decay before Station 4 could cause the event to be misidentified or could cause the reconstructed invariant 
mass to drop out of the D° mass window. About 20% of D° — ► Kir events were lost due to kaon decay. 

D. Dielectron Efficiency 

The trigger efficiency for the decay D° — v e + e~ relative to D° — ► Kir was dominated by the E trigger component. 



As discussed in section [VB, E was used for finding dielectron events while e was used for single-electron events. The 
same prescaled-TFI sample used for the dihadron-trigger-efficiency study was used to determine the efficiency of the 
E trigger. The energy-sum signal of layers E2 and E3, E sum , was digitized by an ADC for off-line study. An efficiency 
curve as a function of the E sum -ADC count was determined and then included in the Monte Carlo. 

Data events with an EM fraction greater than 0.95 and with tracks isolated in the calorimeter were used to 
determine the energy deposited in E2 and E3 as a function of momentum. When an electron was generated from the 
D° — > e + e~ decay, the energy that the electron deposited in the E2 and E3 calorimeter layers was generated based on 
the momentum-dependent (E2+E3)/p distribution. Once the E2+E3 energy stored was determined, an ADC count 
was generated according to the E sum -to-ADC curve. Figure [f(J shows the dielectron efficiency as a function of the 
Monte-Carlo generated D° momentum. The final trigger efficiency for D° — > e + e~ events that passed the Monte 
Carlo geometric cuts was determined to be 60% for both 900A and 1000A runs. 

E. fie Efficiency 

To find the efficiency for D° — ► /i ± e =F relative to D° — > Kir, we used the techniques and tools developed for the 
D° —* fi + 'fj,~ and D° — ► e+e~ efficiency analyses. To adapt to a single-muon event, we used a simple 2-hodoscope 
requirement for the trigger, and demanded 2 out of 3 proportional-tube planes to have hits. The single electron was 
treat ed in a similar manner as the dielectron, the only difference being that the low-energy threshold (see Section 
IX Dj ) required a different trigger-efficiency-to-ADC mapping. As in the dielectron case, the generated electron was 



required to pass the geometric cuts of the Monte Carlo. Figure O shows the turn-on curve for the single-electron 
trigger as a function of the Monte-Carlo generated D a momentum. The resulting single-electron efficiency was 48%, 
and the single-muon efficiency was 74%, yielding a combined efficiency of 36% for the 900A and 1000A runs. 



X. RESULTS 



A. Normalization 



An event was labeled as a D° — > Kir candidate if it was reconstructed as a dihadron event, satisfied the dihadron 
trigger, and passed the impact-parameter and proper-lifetime requirements that were applied to the dilepton decay. 
There was no mechanism to distinguish kaons from pions. However, by Monte Carlo simulation we determined that 
the invariant-mass distribution of D — > Kir events with incorrect particle assignments was much wider than that 
with the correct assignments, 7.1 times as wide for the 900A data sets and 5.4 times as wide for the 1000A set (see 
Figure 112). For each event, the invariant mass was thus computed once as K~ir + and once as ir~ K + . The invariant- 
mass distribution was then fit to a quadratic polynomial for the background and a double Gaussian in the signal region. 
The standard deviation and normalization were allowed to vary for the narrow Gaussian, but the width and relative 
height of the wide Gaussian were constrained to the Monte Carlo values. This condition ensured that the numbers 
of events under the two Gaussian distributions be identical. In addition, both Gaussians were required to have the 
same mean mass. The fit was performed using PAW |2^]. The covariance matrix for the width and normalization of 
the narrow Gau ssia n was then used to find the absolute error associated with the number of reconstructed D° — ► Kir 
decays. Tables IX, fX], and XI give the mean mass, the mass resolution and the total number of D° — ► Kir events 
without any mass cut for each data set. 



B. Signals 



Each data set listed in Tabl e [V| w as analyzed independently. When the impact-parameter and lifetime-significance 
cuts (see Tables VI, VII and VIII) were applied to the dilepton data sample, as shown in Figures [lg|-|l8|, no event 
was found in the signal region, defined as the interval in dilepton invariant mass within which the signal events were 
counted. This interval was ±1.96<r wide and centered at the mean of the Kir invariant mass of the corresponding 
normalization sample. Since the mass resolution depended on th e final-state particles, a different a was used for each 
dilepton decay as computed by MC and tabulated in Table XII. 



C. Summary of Acceptances and Efficiencies 



Tables Kill through XVI give the various contributions to the accept ance x efficiency for the normalization mode 



and each signal mode, together with their errors as estimated in Section X E below. 



D. Branching Ratios 



For each dilepton mode the 900A-Au, 900A-Au-Dedicated-Dilepton, 900A-Be, and lOOOA-Au data sets were com- 
bined to give the branching ratio 



B{D° -> 1+1-) 






x B(D° -> Kir) 



(3) 



where i runs over the three data sets, Si is the number of counts in the signal reg ion i n the dilepton-invariant-mass 
distribution, e, is the efficie ncy fo r the dilepton mode given in Tables XIV , XV and XVI relative to that for the 
normalization mode in Table XIII, and JV, is the number of observed Kir events in the signal region. 

With no signal observed, the upper limit on the branching ratio was determined using the Monte Carlo method. 
A series of branching ratios were calculated according to Equation (ID) . For each calculation the expected number of 
counts in the signal region, Si, was determined from Poisson statistics. Specifically, S was distributed as 



-SoM 



Ml 



(4) 



with M being the actual number of signal counts seen in the data. The quantities Ni, ej, and B(D° — > Kir) were each 
treated as a Gaussian distribution with the same mean and error used in Equation ([|). A minimum of 10 6 events 
were generated for each calculation. From the distribution of the calculated branching ratios, an upper limit on the 
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branching ratio at the 90% confidence level^j was established. By this method, upper limits on the branching ratios, at 
the 90% confidence level, of 1.56 x 10~ 5 for D Q -> fi+ jT , 8.19 x 10" 6 for D Q -> e+e" and 1.72 x 10~ 5 for D° -> fj^e^ 
were obtained. 

We have also employed the method of Cousins and Feldman [J23| , as advocated by the Particle Data Group p4[ , 
for the case in which no signal and no background are observed. In this approach, the upper limits on the branching 
ratio at the 90% confidence level, corresponding to 2.44 events, are 1.65 x 10~ 5 for D° — > n + ijl~ , 8.69 x 10~ 6 for 
D° — ► e + e~ and 1.82 x 10~ 5 for D° — > /i ± e =F decay. These results are about 10% worse than those found by the 



Monte Carlo technique. Table XVII summarizes the single-event sensitivity of our experiment and the upper limits 
for the D° — > l + l~ decays as determined with the Monte Carlo approach and the Cousins-Feldman method. It should 
be noted that the uncertainties in Ni, e,-, and B(D° — > Kir) are not taken into account in the Cousins-Feldman 
method; hence, we favor the Monte Carlo method for determining upper limits. Our D° — ► e + e~ and D° — > n e^ 
limits are comparable to those set by CLEO J4[, whereas the D° — ► /J + /i~ result is about a factor of four worse than 
those of BEATRICE § and FNAL E771 @. 

E. Systematic Errors 

Various systematic errors could affect the results in the branching-ratio determination. We follow Reference p5fl in 
discussing the effect of systematic errors on an upper-limit calculation. 

The uncertainties associated with the p t and xp distributions used in the MC simulation contribute to the error 
in the relative efficiency. This effect was studied by varying each parameter of the pt and xf parametrizations by 
±lcr. The worst case was the variation in xf, resulting in a shift in the absolute efficiency of < 17%. However, in 
each case the efficiency of the dilcpton mode relative to the normalization mode varied much less, only « 2%. The 
temporal variation of the dimuon, dihadron, dielectron, a nd /xe-trigg er efficiencies was investigated by subdividing the 



samples into independent data sets. As shown in Tables XIII XVI, the variations are small. The impact-parameter 



and lifetime-significance cuts, varied by ±lcr, did not have any significant effect on the relative efficiencies. 

XI. CONCLUSIONS 

Three rare or forbidden decays, D° — ► /i + /J~, D a — ► e + e~, and D° — ► /i ± e =F , have been searched for, and no 
evidence has been found for any of these decays. New upper limits on the branching ratio at the 90% confidence level 
are 1.56 x 10~ 5 for D° — > ^ + ^ , 8.19 x 10~ 6 for D° — > e + e~ and 1.72 x 10~ 5 for D° — > ^ ± e =F decay. For comparison, 
the best published limits on these decays are 4.1 x 10~ 6 for D° — ► ^ + /i _ [||J|], 1.3 x 10^ 5 for D a — > e+e" [||, 
and 1.9 x 10 -5 for D° — > /i ± e T B. Our limits for D° — > e + e _ and D° — * ^ ± e =F are the best to date. These limits, 
however, are still many orders of magnitude from the levels at which these processes might be expected to occur. 
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5 That is, 90% of the generated branching ratios were less than this value. 
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FIG. 1. Plan view of E789 spectrometer. 
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FIG. 2. Elevation view of E789 target region. 
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FIG. 3. Schematic diagram of E789 layout from target to beam dump (not to scale). 
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FIG. 4. Isometric view of E789 electromagnetic and hadronic calorimeters. 
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FIG. 5. Conceptual sketch of hodoscope trigger matrix. 
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FIG. 6. Distance between the downstream track and the associated SSD hit for the Y planes in the lower SSD arm for events 
with only one SSD track in either arm. 
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FIG. 7. y- and a;-angle matching between downstream track and SSD track for both SSD arms for events with only one SSD 
track in either arm. 
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FIG. 8. Monte-Carlo generated (left) and accepted (right) distributions of pt, xf and z component of momentum in the 
laboratory frame for D° with SM12 set at 900 A. 
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FIG. 9. Trigger efficiency for dihadron events as a function of pair momentum. The solid curve is the parametrization used 
in the Monte Carlo program. 
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FIG. 10. Efficiency of dielectron trigger as a function of Monte-Carlo generated D momentum. This efficiency was used for 
both 900A and 1000A runs. 
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FIG. 11. Efficiency of single-electron trigger as a function of Monte-Carlo generated D° momentum. This efficiency was used 
for analyzing the 900A and 1000A data. 
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FIG. 12. Reconstructed Kir invariant-mass distribution of Monte-Carlo D° —+ Kit events at 900 A, with entries both for 
correct and for incorrect particle assignments. 
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FIG. 13. Invariant mass distributions for D° — > pi + /i and its associated D° — > if7T distribution for the 900A-Au and dedi- 
cated-dilepton data sets. The cross-hatched area marks the signal region. 



19 





.825 1.85 1.875 1.9 1.925 
GeV/c 2 



FIG. 14. Invariant mass distribution for D° — > /i + fi and its associated D° — > Kir distribution for the lOOOA-Au and 
900A-Be data sets. The cross-hatched area marks the signal region. 
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Kn distribution for the 900A-Au and dedi- 



21 



> 

o 
O 

o 
q 

d 



o 

E 



- 

3330 - 

2960 - 

2590 

2220 - 

1850 - 

1480 - 

1110 - 

740 - 

370 - 

^~ 



> 

u 

O 

q 
d 



1 - 



o 




325 1.85 1.875 



1.9 1.925 
GeV/c 2 




4995 
4440 
3885 
3330 
2775 
2220 
1665 
1110 

50 





.825 1.85 1.875 1.9 1.925 
GeV/c 2 



FIG. 16. Invariant-mass distribution for D —* e e and its associated D 
data sets. The cross-hatched area marks the signal region. 
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FIG. 17. Invariant mass distribution for D° —* /( e T and its associated D° 
cated-dilepton data sets. The cross-hatched area marks the signal regions. 
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TABLE I. Dimensions of targets. 



900A 



1000A 



Length along z (mm) 
Height along y (/jm) 



1.8 

160 



0.8 
110 



TABLE II. Configuration of Silicon Vertex Detector. 



Plane 
No. 


Name 


z Position 
(cm) 


y Position 

(cm) 


View 


Arm 


Number 
of Strips 


1 


Y1B 


-294.54 


-2.125 


Y 


Lower 


316 


2 


Y1T 


-291.36 


0.949 


Y 


Upper 


316 


3 


U2B 


-286.92 


-2.300 


U 


Lower 


372 


4 


U2T 


-283.74 


1.066 


U 


Upper 


372 


5 


Y3B 


-279.30 


-2.758 


Y 


Lower 


436 


6 


Y3T 


-276.12 


1.548 


Y 


Upper 


436 


7 


V4B 


-271.68 


-2.865 


V 


Lower 


500 


8 


V4T 


-268.50 


1.721 


V 


Upper 


500 


9 


Y5B 


-264.07 


-3.364 


Y 


Lower 


572 


10 


Y5T 


-260.88 


2.217 


Y 


Upper 


572 


11 


U6B 


-256.44 


-3.566 


U 


Lower 


628 


12 


U6T 


-253.26 


2.289 


U 


Upper 


628 


13 


Y7B 


-248.82 


-4.018 


Y 


Lower 


692 


14 


Y7T 


-245.64 


2.805 


Y 


Upper 


692 


15 


V8B 


-241.20 


-4.154 


V 


Lower 


756 


16 


V8T 


-238.02 


2.925 


V 


Upper 


756 



TABLE III. Components of TGO Trigger. 



Trigger Name 


Description 


h+h- 


M v ■ M D ■ Su ■ S D ■ NX1 ■ NX3 ■ H 


H + fj,~ 


Mu ■ M D - Su ■ So- NX1 ■ NX3 ■ 2HX4 ■ 2HY4 


e+e- 


Mu ■ M D ■ Su ■ So ■ NX1 ■ NX3 ■ E 


e=F M ± 


Mu ■ M D - Su ■ S D ■ NX1 ■ NX3 ■ e ■ HX4 ■ HY4 


h^e ± 


Mu ■ Mb ■ Su ■ S D ■ NX1 ■ NX3 ■ h ■ e 


h+n* 


Mu ■ M D - Su ■ Sd ■ NX1 ■ NX3 ■ h ■ HX4 ■ HY4 


h ± h ± 


M LIKE ■ (Su + S D )-NX1 ■ NX3 ■ H 


ff 


Mlike -(Su + Sb)- NX1 ■ NX3 ■ 2HX4 ■ 2HY4 


e ± e ± 


M LIKE ■ (Su + S D )-NX1 ■ NX3 ■ E 


e*M* 


Mlike ■ (Su + S D ) ■ NX1 ■ NX3 ■ e ■ HX4 ■ HY4 


h ± e ± 


Mlike ■ (Su + Sd) ■ NX1 ■ NX3 -h-e 


h ± f , ± 


Mlike ■ (Su + S D ) ■ NX1 ■ NX3 ■ h ■ HX4 ■ HY4 





TABLE IV. Average number 


of protons 


i on target and 


triggers 


per 


23-sec spill. 




Run 


Protons on Target 




TFI 






TGO 


TAP 


900A-Be 


3.6 x 10 iu 




2.8 x 10 b 






1.7 x 10 b 


8.7 x W 


900A-Au 


3.3 x 10 1() 




8.3 x 10 6 






1.0 x 10 5 


2.2 x 10 4 


lOOOA-Au 


2.6 x 10 10 




1.7 x 10 6 






1.0 x 10 5 


2.6 x 10 4 
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TABLE V. Summary of data sets. 



Data set 



Protons on target 



AMON ■ SB 



7.2 x lO 13 " 
1.2 x 10 13 
6.6 x 10 13 
2.2 x 10 14 



3.9 x 10" 
8.4 x 10 5 
3.4 x 10 6 
9.8 x 10 5 



TAPS 



3.8 x 10 b 

8.7 x 10 7 
7.0 x 10 7 
7.3 x 10 7 



lOOOA-Au 

900A-Au 

900A-Au-Dedicated-Dilepton 

900A-Be 



TABLE VI. Lifetime Significance and impact parameter cuts for D° - 


-» n+n . 




900A-Au 900A-Au-Dcdicatcd-Dilcpton 


900A-Be 


lOOOA-Au 


Lifetime Significance 1.0 0.60 
Impact Parameter (/an) 25.4 109 


0.90 

78.7 


1.4 
95.3 



TABLE VII. Lifetime Significance and impact parameter cuts for D° 


— » e + e . 




900A-Au 900A-Au-Dcdicated-Dilepton 


900A-Be 


lOOOA-Au 


Lifetime Significance 0.80 0.90 
Impact Parameter (/jm) 48.3 82.6 


0.60 
29.2 


0.80 
55.9 



TABLE VIII. Lifetime Significance and impact parameter cuts for D° 


-> n+e*. 




900A-Au 900A-Au-Dcdicatcd-Dilcpton 


900A-Be 


lOOOA-Au 


Lifetime Significance 1.5 1.3 
Impact Parameter (/im) 190. 97.8 


1.0 

90.2 


1.4 
55.9 



TABLE IX. 


Parameters of D° - 


■+ Kn 


normalization 


for D° - 


■+ (x n search. 




Data set 


900A-Au 




900A-Be 




1000 A- Au 


Dilepton 


Mean Mass (GeV/c 2 ) 

Mass Resolution, a (MeV/c 2 ) 

Total # of D° -> Kn decays 


1.865 

7.0 

606 ± 56 




1.863 

5.6 

1161 ±110 




1.867 

7.3 

971 ± 57 


1.866 

7.9 

2794 ± 396 



TABLE X. 


Parameters of D° - 


-> Kn normalization 


for D° - 


•+ e + e search. 




Data set 


900A-Au 


900A-Bc 




lOOOA-Au 


Dilepton 


Mean Mass (GeV/c 2 ) 

Mass Resolution, a (MeV/c 2 ) 

Total # of D° -► Kn decays 


1.865 
6.1 

707 ± 91 


1.863 

5.9 

1563 ± 154 




1.866 

6.5 

1749 ± 129 


1.865 

6.2 

2867 ± 466 
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TABLE XI. 


Parameters of D° 


-*■ Kn 


normalization for D° - 


■+ fi e T search. 




Data set 


900A-Au 




900A-Be 


lOOOA-Au 


Dilcpton 


Mean Mass (GeV/c^) 

Mass Resolution, a (MeV/c 2 ) 

Total # of D° -> Kn decays 


1.865 

6.0 

255 ±21 




1.864 

6.5 

967 ± 62 


1.867 

7.1 

1429 ± 84 


1.865 

6.3 

2400 ± 337 



TABLE XII. Width of invariant-mass distribution for reconstructed dilepton events relative to that of the normalization 



signal, as determined by Monte Carlo. 






D° -» u+ l x- 


D° -> e+e- 


D° -» M ± e T 


900 A 1.08 
1000A 1.09 


1.09 
1.08 


1.08 
1.07 



TABLE XIII. Efficiencies for D° -> K-k decay. 



900 A 

(3.55 ±0.02) x 10 _a 
0.78 ±0.01 
0.55 ±0.01 

(7.12 ±0.04) x 10~ 2 



1000A 

(2.37 ±0.01) x 10~ 3 
0.82 ±0.01 
0.58 ±0.01 

(9.35 ±0.05) x 10~ 2 



Geometric 

K decays 

Trigger 

Final cuts 



Total 



(1.08 ±0.02) x 10" 



(1.05 ±0.02) x 10" 



TABLE XIV. Efficiencies for D° ->■ fi + ^ decay. 





900A 


1000A 


Geometric 

Trigger • ID • isolation 

Final cuts 


(3.97 ±0.02) x 10" d 

0.36 ±0.01 
(7.30 ±0.04) x 10~ 2 


(3.05 ±0.02) x 10 -3 

0.50 ±0.01 
(9.85 ±0.05) x 10" 2 


Total 


(1.05 ±0.02) x 10~ 4 


(1.51 ±0.02) x 10" 4 





TABLE XV. Efficiencies for D° -► e+e" 


decay. 






900A 




1000A 


Geometric 

Trigger • ID 

Final cuts 


(3.87 ±0.02) x 10~' 5 

0.60 ± 0.01 
(7.15 ±0.04) x 10~ 2 




(3.28 ±0.02) x 10" J 

0.60 ±0.01 
(10.0 ±0.05) x 10" 2 


Total 


(1.65 ±0.03) x 10~ 4 




(1.98 ±0.04) x 10" 4 
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TABLE XVI. Efficiencies for D° 



±oT 



ix e 



decay. 





900A 


1000A 


Geometric 

Trigger • ID 

Final cuts 


(4.18 ±0.02) x 10 -3 

0.36 ±0.01 
(7.04 ±0.04) x 10~ 2 


(3.42 ±0.02) x lO" 3 

0.36 ±0.01 
(10.0 ±0.05) x 10" 2 


Total 


(1.06 ±0.03) x 10~ 4 


(1.28 ±0.04) x 10 -4 







TABLE XVII. Summary 


of results. 




Decay 
Mode 


EiNid 


Single-event 
Sensitivity 


Limit from MC 
calculation 


Cousins- Feldman 
Method 


D" - e+e~ 
D° - M ±eT 


5830 
11100 
5300 


7.34 xl0 -b 
3.84 xl0~ 6 
8.08 xl0~ 6 


1.56 xl0~ b 
0.82 xl0~ 5 
1.72 xl0~ 5 


1.65 xl0~ b 
0.87 xl0~ 5 
1.82 xl0~ 5 
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